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ABSTRACT 

The effects of a neutron charge distribution are considered. 
It is suggested that the charge distribution is of the form 

i 

p-TY^ ^ 

The effects of this charge distribution for various values 

of the parameters c and f , on the mean square radius of carbon, 

and on the charge form factor of carbon at high q are shown. The 

effect of the neutron becomes appreciable in carbon above the 

2 

diffraction minimum. The neutron in this range of q adds approx- 
imately 50% to the carbon form factor. 
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1. INTRODUCTION 



A review of current literature concerning the nuclear charge 
distribution shows that the effect of the neutron is largely ignored 
except in deuteron studies. This paper proposes to suggest what 
effect, if any, a finite neutron will have on the distribution of 
charge in any nucleus. 

Before the charge distribution in a nucleus can be calculated, 
the wave functions describing the motion of the nucleons and the 
nucleons' internal charge distributions must be known. The proton 
charge and magnetic moment distributions and the neutron magnetic 
moment distribution are known. The neutron charge distribution is 
not well known, however . 

If wave functions are known, then the neutron charge distribution 

may be related to the measured nuclear charge distribution. In this 

paper the diffraction minimum observed in electron scattering by a 
12 

C is used to determine the nucleon wave functions. 
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2. THE NUCLEAR FORM FACTOR 



In order to develop a theory for the charge distribution 
of a nucleus, it is assumed that each nucleon maintains its 
identity in the nucleus. This assurrption leads to models which 
successfully predict various experimental results. The assurrption 
iirplies that each nucleon in the nucleus may be described independ- 
ently of the other nucleons. The nuclear wave functions may be 
written as follows: 



( 2 - 1 ) 




) 

A l 



y'it) vW) 






r^ describes the position of the center of mass of the 
th 

i particle. The ^ 1 inplicitly include a spatial wave function, 

hereafter referred to as 'V , spin and isospin functions . All of 

these functions are orthonormal. The use of the othonormality 

will be implicit in what follows. The total wave function, $ , 

is also normalized and is antisymmetric on interchange of any two 

nucleons (a consequence of the nucleons being fermions) . 

The purpose of this paper is to examine the nuclear charge 

distribution p (R) . This cannot be measured directly. What can 

be measured is the differential cross section for scattering electrons, 
d<T 

asi' 
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The cross section may be separated into a term due to the 

charge distribution, and a term due to the magnetic moment of the 

nucleus [1] . The charge distribution term leads to the charge form 
2 

factor, F(q ), which, m Bom approximation, is the Fourier trans- 
form of the normalized charge distribution. 

<2-2) c/ 3 r 

The Bom approximation will be used throughout this paper. 

The nuclear charge distribution (R) , is determined by the 
nuclear wave functions and the charge distribution of the nucleus. 

(2-3) s ( /.'lei 3 ? 

-A 

where is the nucleon charge operator: 

<2 * 4) / -- 

The f^ are the charge operators for the individual nucleons. 
Putting equations (2-3) and (2-4) into equation (2-2) yields 



(2-5) 



Kf) = 






Let <|^ (r . ) be denoted by ^ . Note that after integrating over all 



except the i 1 coordinate, we have 



( 2 - 6) /£*{(*- 



T. /3A 

r - 



A 



■ 
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We see from equation (2-5) and (2-6) that, if the single 
particle wave functions are orthogonal, the charge distributions 
for all of the protons are equal, those of all the neutrons are 
equal and the charge distributions are additive. 

Using (2-6) in (2-5) implies 



<2-7) F(f) 

C c = i n J 5 ' J 



The relation between R and r. is R = r. + r . ' • 

l ii 

Insert the following equation into (2-7) : 




Equation (2-7) becomes 







Integrate over R, 



( 2 - 10 ) 



R?0 - 



J_ £ -L 





Separate integrals into one over A and one over r\ ' : 
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Split the sum into one over protons, f^ = f , and one over neutrons, 



f. = f 
1 n 






Define /F = 







( 2 - 12) 1 



p = 

F *' 



<*■*> (j= . 



Since we have integrated over r, F^_, = F^ for all i. Use (2-12) 
and (2-13) in equation (2-11) . 



(2-14) 



no’twns'&n! 



Define 
(2-15) P 



£ K 



(m-A* A j - 1 ^ 



Then write equation (2-14) as 
(2-16) 



Rf) = [ f> + i ^ • 



( 1 ) 



The nuclear form factor has been separated into three parts: 

F , the form factor of the proton 
P 
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(2) F , the form factor of the neutron, and 

n 

(3) a form factor due to the distribution of the centers of 
mass of the nucleons. 

F'cm-A is valid for any nucleon distribution,^ (R) ; however,^ (R) 
may be expanded (since any function of angle may be expanded in terms 
of spherical harmonics ) as follows: 



charge distribution. This yields an electric monopole term plus higher 
order multipole terms in the form factor. The simplest nuclei to 
analyze are then the spherically synmetric nuclei, since these have 
only the electric monopole term. All further work in this paper will 
deal only with this type of nucleus. 



(2-18) /?($ ) - 




where ^ is f* (R) averaged over all angles or the spherically 



symmetric part of 
2 

F(q ) may be expanded in terms of the multipole expansion of the 
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3. MEAN SQUARE RADIUS DETERMINATION 



A quantity which is experimentally easy to determine is the 
mean square radius. We show below that, theoretically, this radius 
may be expressed in terms of the particle wave functions and in terms 
of the charge distributions of these particles. 

The form factor for a spherically synmetric charge distribution 
is _ 

(3-D jV v /^ r WV 

. iJ Jclr . 

o / 



We have on expanding sin (qr) : 




This implies that 




From equation (2-16) we have 



(2-16) F(f) = [ f 0 + -e" ^ 
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Using this in equation (3-3) yields 



(3-4) <r"> 





i d? 









J r J 



V“ 




In equation (3-4) we have used 




The directly measurable quantities in equation (3-4) are 
the mean square radii of the nucleus and of the proton. The mean 
square center of mass radius is dependent on A and the configuration 
of the nucleons in the nucleus. This configuration may be different 
for different elements of equal A. 

It is thus possible to write the nuclear form factor in terms 
of the nucleon form factors and nucleon wave functions which describe 
the states of the nucleons. See equation (2-16) . However, starting 
with experimental data on the nuclear form factor, it is impossible 
to separate the contributions due to the nucleon charge distributions 
and that due to the nucleon motion unless a model for the nuclear or 
nucleon structure is assumed. That is, the analysis will be model 
dependent. The validity of any model is measured by how well it 
agrees with experiment. 
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4. POSSIBLE NUCLEON CHARGE DISTRIBUTIONS 



There has been a great deal of work done on the proton 
which suggests that it has a gaussian charge distribution [2] . 



(4-1) 



to • 



-*• 



ci 



The form factor corresponding to this charge distribution is 

i>JL 

(«> W r £ ' ** 



Figure 1 shows a corrparison of this form factor to experimental 



data [2] . We assume here that f is correct, and that the charge 

P 

distribution and the matter distribution are identical. With this 
assumption it seems valid that the matter distribution in the neutron 
could also be distributed in the same manner. That is to say, the 
matter distribution of the neutron is a gaussian with approximately 
the same radius as the proton. It is known that the total charge 
of the neutron is zero. However, its magnetic moment suggests that 
the neutron does have a charge distribution. By analogy with the 
proton the following form is appealing: 



(4-3) 



- --h 

/ I fi _J ^ -f u 

r ~ t 34 



This is a double gaussian, each with unit charge. The negative charge 
has the larger radius. 

Our assumption that the matter distribution of the proton is the 
same as that of the neutron implies that the radius of the neutron's 
negative (larger) gaussian is approximately that of the proton. 
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V/h 




Figure 1 
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We take f = .44F (slightly larger than the proton) . This leaves 
one parameter to fix in the neutron. The free parameter will be 
fixed in this paper by using the harmonoc oscillator model and 
experimental carbon data. (For alternative method see Appendix C.) 

The form factor correspond to f R is, in Bom approximation 




, 2 

Figure 2 is F n for various values of c compared to experimental 

data [3]. This data does not deny any of the suggested curves, but 

favors F =0. 
n 
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5 . HARMONIC OSCILLATOR WAVE FUNCTIONS 
2 

In order to determine < r 7 we are forced to assume a model 

n 

2 

to predict ^ r^'’ . The harmonic oscillator potential is a simple 
potential which appears to agree well with experiment. This theory 
has the added advantage that its parameters are fixed independently 
of the mean square radius of the charge distribution. 

Consider non-interacting particles in a spherically syrrmetric 
harmonic oscillator potential well: 

( 5 - 1 ) V(. r ) - /auJ X Y X ‘ . 



The stationary states or orbitals occupied by the particles are 
solutions of Schrodinger 1 s equation. 

1 



(5-2) 



[. - e ^ 



4 ^ 



Using equation (5-1) , and defining a = 2 / equation (5-2) can be 

m w 

written as 



(5-3) 



I ± 



1 r 



A 

z. ^ 



_ Z ”W(?) = E ) 

^ ^ J 



(5-4) where 
i I 






^ 5 7'Tr ^ 3 



+ , V I v 



Assume 



(5-5) 



l^(_ r) - VsC?) y . 



Using equations (5-4) and (5-5) in equation (5-3) 
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(5-6) 



r or 



or ^ or i 





i _ 

S 



0 ($ *)£>) 







j'(j + 0 

The solutions to the angular equation are the spherical harmonics. 
(See Appendix B.) 

For convenience in solving the radial equation, define 

6 Hr) 

Y“~ 



(5-7) *R(r) 



(5-8) 



The radial portion of equation (5-6) may be written as 

^ | lil) / b - y wto 1 f 1 £f v( *')V ] L /(k) - o 

A r‘ m‘ L / 



The solution of equation (5-8) and equation (5-7) yields: (See 



Appendix B.) 



* 



1 “ 6 C< Ch" 



(5-9) D = Y c 
lx A»< 



< -0 



lr 



i < 



These equations also yield: 

(5-10) = ^to [in t .(•' + 1. ] 

Where A = 2N + 1 denotes the energy levels. The harmonic wave 
functions are then (See Appendix B.) 

(5 - n) to,.'- . 
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The above solutions are based on the central potential having 
a fixed origin. In reality, the central potential should be measured 



from the mean position of the nucleons. It has been shown that this 
corrplication introduces a factor <Z * IA into the form factor [4]. 

From equation (2-4) and the above consideration we have (See 
Appendix B.) 




From equation (2-16) we have 



n, 1, m must be known for each nucleon) are known, theh the total 
nuclear form factor may be written. For exanple, for p-shell 

nuclei (4CA^16, 2^ 8) [5] : 



The F , from Appendix B have been used. In Carbon A = 12, Z = N = 



(5-13) 




If the harmonic oscillator model is applicable and if the 

nucleon form factors F and F and the nuclear configuration (i*e., 

n p 
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to|> 



Or the total form factor is 



%vt) * (*> + f- 



no- *T) 






'O’O 



2 cr 

F has a zero at q = 2_ . Therefore, our cross section measurements 

aT 

should also shew a zero at this point. Experimental data shew a 

2 

marked minimum in the cross section which can be used to fix a . 



6. COMPARISON TO EXPERIMENTAL DATA 



In the harmonic oscillator model, the center of mass form 



factor, F cm _ A / is completely determined by the position of the 

2 2 2 
minimum in a plot of F (q ) vs. q . 

Carbon offers a very good method of fixing the neutron parameters 
since it is spherically symmetric and the experimental data show a 
prominent diffraction minimum [6] . 






( 6 - 1 ) 



( 6 - 2 ) 



A ' - JL -- f-, F -> - 2 

<■* c l 3.11 r 



Figure 3 is a comparison of equation (6-1) to experimental values 
for the neutron parameters of Figure 2 [6] . This figure indicates 

that the best neutron form factor is given by 

/< . - . .\? 



(6-3) 



- . 
F: £ - C 

1 Y\ 



l 



(Jj ) 1 



A i 



The mean square radius of carbon is then 



1 ^ \ 
< r > ( ; 



(6-4) 



. ,(,* F* - ,u F’ + 5,41 F'.- (,.01 F 1 

This radius is, of course, model- dependent because of the manner 
2 

in which <r has been found. However, it agrees well with current 

an-A 3 



model independent determinations 



of<r 2 ^( 



c at Darmstadt in which 



^ r 2 ^ c = 



.587 + .19 F [7] . 
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7. DISCUSSION 



The preceeding results suggest two alternative approaches 
for determining the validity of the neutron radius in equation 
(6-4) . First, assume that the nucleons are described by harmonic 
oscillator wave functions, in which the range parameter a is adjusted 
to give the experimentally observed diffraction minimum in F (q ) . 

This allows us to write 



(7-1) 



iL < Y-; ? r < r* > 

■c 



- < 



- 



< r 1 > 



We would then use the model independent determined at low q , 

2 1 

the accepted value of r^ determined from H x data, and the shell 
model prediction of < r X . Extreme accuracy in the measured 

du — A 

2 

quantities is required due to the small size of < r n > . 

The second alternative would be to assume that the nucleon 
center of mass distribution is the same for elements with equal A 
but different Z. This assumption is based on the charge independence 
of nuclear forces. The assumption could be checked in a model- 
dependent way by seeing whether or not their diffraction minima 
occurred at the same value of q. Measuring the model- independent 
radii of two elements with the same atomic weight at low q and 



subtracting the following: 



-ch*. r (f- . g' => 



(7-3) 









A (.*-£') 



<< r 
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Thus^-r^ is obtained in terms of only the model independent 

2 2 ^ 

mean square radii. The assunption, however, that^ r ' A r 



om-A 



would only be valid if the nucleon configurations of the elements 

were identical. One possibility would be to use chramium-50 and 

titanium-50. These two nuclei have total angular momentum (J) equal 

to zero. Therefore, they have only a monopole form factor. Any 

other pair of mirror nuclei would do equally well. 

Both of these methods could be checked by comparing the 

theoretical form factors to the experimental form factors at high 
2 

q . This is the manner that the neutron parameters were fixed m 
this paper. 

The limitation on the parameter f in the neutron charge 
distribution (that the neutron matter distribution equals the 
negative charge distribution) is of questionable validity. If 
this restriction is removed, and we use instead the electron- 
neutron interaction condition. 



(7-5) 



0 . \ 

J 






which implies 



(7-6) < r * } = - , iz& F 1 



2 

Then the difference between the two neutron parameters is .084 F . 



(7-7) 
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F.« £ 



When this restriction is used and the form factor 



is again fitted to various 



experimental data, a good fit is obtained for c = .09 F , 

2 2 
f = .174 F z . 



This choice of parameters leads to a form factor which gives 

2 2 -2 

a smaller radius at low q , but at q of about 25 F the radius 

is much greater than with the previous choice of parameters. The 

form factor of the neutron, with this choice of parameters, becomes 

2 -2 

about three times that of the proton at q = 25 F . This doesn't 
seem reasonable. The previous selection of neutron parameters is, 
consequently, preferable. 
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8 . CONCLUSION 



Current experimental data point to the possibility that the 

charge distribution of the neutron has a greater mean square radius 

2 2 

than previously thought. It appears that a value of< r n ^ = -.2F 
is not unlikely. In fact, such a value is supported by current 
measurements of nuclear radii. 

It should be noted that it is presently inpossible to separate 
the effect of the neutron from the effect of the wave functions in 
currently available experimental data on electron scattering. A 
modification of the wave functions can give form factors which fit 
the data over a large range of q without considering any neutron 
effect. 

The approach suggested in this paper has the advantage of 

simplicity. Simple wave functions and a simple neutron give agree- 
2 

ment with F(q ) data. 

If the charge distribution of the neutron can be determined 
independently (for example, through neutron-electron interactions) , 
then this information could be used in fitting wave functions to 
experimental data. 
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APPENDIX A. COMMON INTEGRALS 
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APPENDIX B 



HARMONIC OSCILLATOR WAVE FUNCTIONS AND FORM FACTORS 



1 . Spherical Harmonics 
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2 . Radial Wave Functions in n j. 
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3 . Harmonic Oscillator Wave Functions 

T _ i_ Y t 



d/ = ^ 7 

Tu 0l ^ ' 



(( ' -2 fVv 



K. - [ 



0^ Ti Vt 



t 7-Tw 



\ £o$6 d ^ 



‘f'o.o r [ 

( r 1 I'f i ^ --£\ 

To.ii - L <S r V* ~ J ^ J<AJ6 c \ 

,, _ r 3. H . 2 -£v 

'»•» " L 2. TT Vt J C » - -i ^ ^ 

[ aV^rC"]’ ^‘d' A ‘ (3^'.i) 



+ 



<?e i \ 



n-. 



~ [ 
lotjri L 



«’-n* j 






0- _. 

r l e * n <^>5 6? J /^6) ^ 



ic<f 



2 <5 Tt 



e 



v* 1 



JV^ £> C 



12^ 



4. Electric Monopole Form Factors for Independent Particles 
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APPENDIX C 



CLASSICAL MAGNETIC MOMENT DETERMINATION 



Classically the magnetic moment m is defined as 



(C-l) |V\ = 



\ 






j' r X 0 C r)c\ ^ r 



[ 10 ] 



where T is the distance from sore point of interest and J is the 
current density. For an arbitrary charge distribution rotating 

A 

about an axis with velocity < 

(C-2) ^ -s (. Uy> K V “" ) C r ' - 

Using this in equation (C-l) yields 

L — • ' ' (r) a 



(C-3) - 



3 



it' r * UKr r 



Ir* 



Using the vector identity r X X T r (f rUo - ( r ^ ) Y' 
this becomes 



(04) 






= ~ic n r!T5 “ v ‘' w ^ cl 



Vie are interested in finding the magnetic moment for neutrons 
and protons projected onto the R axis. For a spherically symmetric 



charge distribution we have 



(05) 



fA • YC -=■ 

£ C 



jV 1 tv) cl >' - c os cfr 



The second term integrates to zero when the angular integration is 
performed. 



So the magnetic moment about the Z axis is for the proton 



(C-7) M 






L^y r 



£C 

and for the neutron 



(C-8) |\A 



? y i" 



u) ^ 

rz 



T S 
1 P 






The ratio of the proton moment to the neutron moment is then 



(C-9) K 



yV\^ f 



1,7-1 

l.*\\ 







assuming £*) - tJ) f> 



This agrees in sign with the value predicted fori r^ /" as 
2 

coirpared to i r p > but is only the correct order of magnitude and, 
therefore, gives neutron parameters larger than those in this 
paper. 

(C-10 ) < *1 7? “ 
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at high q are shown. The effect of the neutron becomes appreciable in carbon above 

2 

the diffraction minimum. The neutron in this range of q adds approximately 50% 
to the carbon form factor. 
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